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Abstract 

Steady forcing at the wall of a channel flow is studied via DNS to assess its ability of yielding 
reductions of turbulent friction drag. The wall forcing consists of a stationary distribution of 
spanwise velocity that alternates in the streamwise direction. The idea behind the forcing builds 
upon the existing technique of the spanwise wall oscillation, and exploits the convective nature of 
the flow to achieve an unsteady interaction with turbulence. 

The analysis takes advantage of the equivalent laminar flow, that is solved analytically to show 
that the energetic cost of the forcing is unaffected by turbulence. In a turbulent flow, the alternate 
forcing is found to behave similarly to the oscillating wall; in particular an optimal wavelength is 
found that yields a maximal reduction of turbulent drag. The energetic performance is significantly 
improved, with more than 50% of maximum friction saving at large intensities of the forcing, and 
a net energetic saving of 23% for smaller intensities. 

Such a steady, wall-based forcing may pave the way to passively interacting with the turbulent 
flow to achieve drag reduction through a suitable distribution of roughness, designed to excite a 
selected streamwise wavelength. 
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I. INTRODUCTION 



In recent years several attempts at controlling turbulence through a number of wall-based 
forcing methods have been reported often aimed at frictional drag reduction (DR). A 
large number of works, exploiting both numerical and experimental approaches, has been 
devoted to this goal. 

In this paper we focus on spanwise wall-based forcing, i.e. on a class of forcing methods 
designed to modify favorably the turbulent flow by introducing an external action directed 
in the spanwise direction. Early work addressing the modification of wall turbulence by 



creating a cross flow can be traced back to Bradshaw & Pontikos Spanwise forcing of 
turbulent flows has been reviewed by Karniadakis & Choi [3]. In 1992, Jung, Mangiavacchi 
& Akhavan [4j introduced the spanwise-oscillating wall technique, where the wall of a fully 
turbulent channel flow is subject to an alternate harmonic motion in the spanwise direction. 
If W indicates the spanwise velocity component at the wall, the law that defines this forcing 
method is: 

W = Asm(^t] (1) 



x T 

where t is time, and A and T are the oscillation amplitude and period, respectively. Once A 
and T are set within the optimum range from the viewpoint of DR performance, the authors 
observed, by means of direct numerical simulations (DNS), that a strong suppression of 
turbulence occurs in the wall region, accompanied by a significant reduction of the mean 
friction. 



The analysis has been carried on in successive studies: we recall Quadrio & Sibilla 
for the pipe flow, and J. -I. Choi, Xu & Sung and Quadrio & Ricco [7, sj], as DNS-based 
analyses which contributed new datasets and detailed descriptions of the flow. Laboratory 
experiments, due to Laadhari et al. [9j and K. S. Choi id Ill], complemented the numerical 
works and addressed the issue of dependency of DR on the value of the Reynolds number. 
Alltogether these works have contributed to demonstrating that the natural friction drag 
of the turbulent flow can be reduced (at least at moderate values of Re) up to 45% for 
A + m 25 (quantities with the + superscript are made dimensionless with viscous wall units). 
An optimal value of the oscillation period exists, namely T^ t = 100 — 125, that yields the 
maximum DR at all amplitudes. The net energy saving, that substracts from the reduced 
flow-driving pumping power the power expense required to move the wall against viscous 
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resistance, has been addressed first by Baron and Quadrio [12| . and today it is recognized 
that it can reach up to 10%. Crediting a commonly accepted qualitative explanation, the 
transverse oscillating boundary layer induced by the wall motion explains the drag reduction, 
since it produces a phase displacement between the wall-layer turbulence structures, capable 
of weakening the viscous wall cycle. When phase-averaged, this time-alternating layer in 
the turbulent regime has been found to coincide with the oscillating laminar Stokes layer, 
: 'or which an analytic solution is known as a classic solution of the boundary layer equations 
l^ . This has been recently exploited to determine a parameter 8] that is capable of scaling 
linearly with DR and thus makes it possible to predict DR capabilities. 

Advantages and drawbacks of the oscillating-wall technique are obvious. It presents 
energetic performance that could make it worth of practical implementations, and thanks 
to its open-loop character it does not need distributed sensors or actuators, that would be 
still unpractical with the technology available today. On the other hand, by its very nature 
this technique requires moving parts, and thus does not lend itself to be implemented as a 
passive device, which is on the other hand the most appealing possibility application-wise. 

In this paper, we aim at extending the oscillating-wall technique, to take one further 
step towards the long-term objective of a successful practical implementation; in particular 
we want to translate the time- dependent forcing law expressed by Eq.([T]) into a stationary 
formula. This goal can be achieved by exploiting the convective nature of wall-bounded 
flows. Though at the wall the mean velocity profile annihilates, it is well known that the 
convection velocity of a turbulent wall flow (or, more precisely, the convection velocity of 
turbulent fluctuations) resembles the mean velocity profile only in the bulk of the flow. Kim 
& Hussain jljj have shown some years ago that near the wall, say below y + = 15, the 
convection velocity becomes essentially independent upon the wall distance, and remains 
constant at the value ~ 10. It is this near-wall value of the convection velocity that will 
enable us to translate the temporal forcing from Eq.(0Q) into a spatially-oscillating forcing, 
yielding the following forcing law that is expected to modify the turbulent flow in a similar 
manner to the oscillating wall: 

W = Asm(^x). (2) 



Here x denotes the streamwise coordinate, and X x is the forcing wavelength. The result- 
ing distribution of spanwise wall velocity is sketched in figiH where the coordinate system 



W=Asin(k x) 



FIG. 1: Sketch of the wall spanwise forcing discussed in the present work. 

employed in this paper is also indicated. 

Though the control law ([2]) has never been considered in the literature, one paper where 
a similar space-time extension has been discussed in the past is the one by Berger and 



coworkers 15J . In a parametric DNS-based study a spanwise-oriented Lorentz volume force 



was simulated to obtain the following two forcing configurations: 



Be^shx^t), (3) 
Be~ y/A sin fer ) , (4) 



where A is the penetration depth of the forcing, and B its intensity. For certain values of 
the parameters, the streamwise-dependent law (HI) produced a DR comparable to that of the 
time-dependent law ([3]), with the added benefit of an improved energetic balance. However, 
the difference between a body force and a wall-based forcing is substantial, as discussed for 
example by Zhao et al. flfl. and no definite conclusion can thus be drawn a priori from the 
study by Berger et al. [15] with respect to the oscillating velocity, which is strictly a wall- 
based forcing. Indeed, it will be shown later that some 

(for example the existence of an optimal wavelength that depends on the forcing intensity) 
do not apply at all to the type of forcing considered here. Schoppa & Hussain [l7|] showed 
that a significant amount of drag reduction can be achieved by introducing in a channel 
flow a spanwise velocity gradient, generated in the DNS by a large-scale, streamwise-aligned 



and x-independent rolls. Such velocity gradient weakens the near-wall cycle by suppressing 
the transient growth of streaks that would otherwise be stable according to normal-mode 
analysis. This mechanism is then addressed in a subsequent paper [18j, where the role of 
streamwise- varying spanwise perturbations is highlighted in the context of the turbulence 
regeneration cycle. The dominant wavelength of the streak waviness is found to be very 
similar to the streamwise length of the coherent structures educed from conditional analysis 
of turbulent flow fields [lij, namely about 300-400 viscous units. 



Aim of the present paper is to investigate the DR and energetic performance of the 
steady control law (j2J). Thanks to an accurate data set, purposely obtained by several DNS 
simulations, the X x — A parameter space will be explored in detail. The value of the Reynolds 
number will be fixed at Re T = 200 (based on the friction velocity of the reference flow and 
half the channel gap); any dependence on this flow parameter is not discussed here. The 
effects of the forcing j2} on a laminar channel flow will be preliminarly studied. We will see 
that a laminar solution, although approximated, can be obtained in analytical form, and 
this solution will then be used to help understanding the DR properties of the forcing (J2J) 
when used on a turbulent flow. In the end, we will be able to connect this forcing to a 
physical application that does not necessarily involve a moving wall. 

The paper is organized as follows. SecJTT] contains a theoretical discussion of the associ- 
ated laminar flow, that is useful for predicting some of the energetic characteristics of the 
present technique (further analytical details are deferred to an Appendix). In Sec JIIII the nu- 
merical simulations of the turbulent case and their discretization parameters are described. 
Results about DR and energetic performances are reported and discussed in Sec JIVl together 
with a visual and statistical description of the forced flow, compared to the reference unper- 
turbed one. Some conclusive remarks, including examples of how the present forcing can be 
implemented in practice, are given in the concluding SecJVl 

II. LAMINAR FLOW 

We consider first the laminar flow in a plane channel subject to either the temporal 
boundary forcing (TjQ) or the spatial boundary forcing (j2j). For the temporal case, it is easy to 
show that the incompressible momentum equation for the spanwise component w decouples, 
so that the streamwise flow is described by the classic Poiseuille parabolic solution, and 
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FIG. 2: Contour of the w field for the laminar TSL and SSL. Top: temporal forcing w(0,t) = 
sin(oji). Bottom: spatial forcing w(0,x) = sin^Kx). Contours at (—0.45,0.1,0.45), dashed lines 
indicate negative values. 

the entire flow consists in this parabolic profile plus a spanwise alternating motion, that is 
identical to the oscillating transversal boundary layer that develops in a still fluid bounded 
by a wall subject to harmonic oscillation (i.e. the so-called Stokes second problem, that 
possesses a classic analytical solution |13|). This oscillating boundary layer will be called 
Temporal Stokes Layer (TSL) in the following. 

We address now the laminar flow subject to the space- varying boundary forcing (j2J), to 
define and discuss the spatial equivalent of the TSL, that will be called Spatial Stokes Layer 
(SSL). We leave some analytical details to Appendix IVHl but the main ideas are given here 
as follows. 

In both the TSL and SSL, the field of the spanwise velocity component w is a function of 
two independent variables, namely t, y for TSL and x, y for the SSL. We start by observing 
in figj2] how closely these w fields for TSL and SSL resemble each other. The abscissa in 
figj2]is t/T for the TSL, whereas it is changed to x/X x for the SSL. Both fields are computed 
with the DNS code, although the analytical expression of the former is available [lsI ]: 

w(t,y)=C t ?St[e iut e^ 5 *], (5) 
where C t is a (real) normalization constant, u = 2tt /T is the angular frequency of the 



oscillation and St is the thickness of the TSL, defined as: 




(6) 



Figure [21 where the parameters T and are in the range of interest for turbulent 
DR purposes, illustrates the wall-normal structure of the TSL and the SSL. Within this 
parameter range, both Stokes layers are thin compared to the channel half-width h. The 
convection velocity being approximately constant over such a small wall-normal extension 
explains why the two contours in figure [2] look very similar. 

An analytical solution for the SSL can be arrived at under the small thickness approx- 
imation. After some analytical efforts, described in the Appendix, the w(x,y) field of the 
laminar SSL can be shown to obey an expression, similar to Eq.fjSJ), that contains an Airy 
function instead of an exponential function: 



where C x is a (real) normalization constant, k = 2ir/\ x is the forcing wavenumber, Ai is the 
Airy function of the first kind [20j], and 5 X is the thickness of the SSL, defined as: 



evaluated at the wall. The w(x,y) laminar field computed by DNS and plotted in figfj] 
(bottom) is virtually indistinguishable from the same field as computed from the analytical 
solution (j7j). 

Comparing the two expressions ((5} and ([7]), a qualitative difference between the SSL and 
the classical TSL can be observed: the former is not decoupled from the longitudinal mean 
velocity profile u(y), but it depends on u Vy o through the thickness 5 X given by (jSJ). This 
remains without consequences in the laminar case, where the streamwise flow is not affected 
by the wall forcing, whereas it will become important in the turbulent case. 

A. Comparison with the mean spanwise turbulent flow 

It is well documented j3] for the time-oscillating wall that the phase-averaged (w) xz (y) 
profile is identical to the laminar solution expressed by Eq.(jS]), except for the initial transient 
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(8) 



In this expression, u y fl represents the gradient of the streamwise mean velocity profile 
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FIG. 3: Comparison between laminar and turbulent mean w fields for temporal and spatial forcing. 
Left: analytical solution (lines) of the TSL and turbulent (w) field (symbols) computed by DNS. 
Right: analytical solution (lines) of the SSL and turbulent (w) tz field (symbols) computed by DNS. 
Profiles are taken at different (temporal and spatial) phases during the oscillation. The laminar 
solutions are linear in A, whereas the turbulent profiles are for A + = 12. Turbulent temporal 
case is at Re T = 200 and T + = 125. Turbulent spatial case is at Re T = 200 and A+ = 1250. 
To emphasize collapse with proper scaling, the SSL analytical solution is computed for a different 
wavelength, namely A^T = 1900, yielding a different 6 X , so that the role of y/5 x as similarity variable 
is highlighted. 

where the oscillation is started from rest. (The operator (-) xz indicates averaging along the 
homogeneous directions x and z.) This is well illustrated by figJ3] (left), where the wall- 
normal distribution of w in the TSL after Eq.(jSJ) is plotted at various phases during the 
cycle, and compared with the turbulent (w) field. The agreement, which has been related 



by Ricco & Quadrio [2l| to the vanishing y, z component of the Reynolds stresses tensor, is, 
as expected, excellent. 

The same result is shown in figE] (right) to hold true for the SSL case. Now, of course, the 
analytical solution (j7j) has to be compared with the turbulent (w) tz field. The role oiy/6 t as 
similarity variable is well known for the TSL. To emphasize that y/S x plays the same role for 
the SSL, the comparison between laminar and turbulent profiles is carried out for two cases 
with different k and thus, according to (jSJ), with different 5 X . The agreement between the 
two profiles is excellent. This confirms the analogy between TSL and SSL in terms of the 
connection between the laminar solution and the phase-averaged turbulent solution. This 
analogy can be used (see lIVp to increase our predictive capabilities when the spatial forcing 
(J2J) is applied to a turbulent flow. 

III. NUMERICAL METHOD FOR DNS OF TURBULENT FLOW 

We turn now to the turbulent case, that will be dealt with by using the Direct Numerical 
Simulation (DNS) technique. The computer code used to solve the incompressible Navier- 
Stokes equations is a parallel DNS solver, based on mixed discretization (Fourier expansion 
in the homogeneous directions, and compact fourth-order accurate explicit compact finite 
difference schemes in the wall- normal direction), recently developed by Luchini and Quadrio 
22)]. The turbulent flow in an indefinite plane channel is considered, with gap half- width 
given by h. One reference simulation without forcing has been carried out at Re T = 200, 
where Re T is the friction Reynolds number, defined based on h and the friction velocity u T , 
that will be used as reference velocity throughout the paper unless otherwise noted. The 
size of the computational domain is L x = Qnh, L y = 2h and L z = 3nh. The total averaging 
time is rather large, to allow for well converged time-averaged results, and amounts to 
approximately 10 4 viscous time units. The employed number of modes / points is given 
by N x = 320, N y = 160 and N z = 320, that yields a standard spatial resolution. For 
the reference simulation a value of Cf = 7.94 • 10 -3 (defined based on the bulk velocity) is 
obtained, in agreement with the correlation Cf — 0.0336i?e~ ' 273 given by Pope [231 ] . 

Some 40 simulations with Eq.([2j) used as boundary condition at the channel walls have 
been then carried out for different values of X x and A. Owing to the periodic boundary 
conditions employed in the homogeneous directions, an integer number of forcing wavelength 
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Aa; must be contained in the domain length L x . The systematic variation of X x has thus 
required, mostly for the largest values of X x , slight adjustments of L x . Compared to the 
reference value of L x = 6irh, 3 cases had an actual value of L x = 7nh, and one was at 
L x = 5irh. Since such changes are of limited entity, the number of streamwise Fourier modes 
has been left unchanged, with the implied little change in streamwise spatial resolution 
for these few cases. Properly accounting for the initial transient, where the wall friction 
decreases from the unperturbed value to the reduced value given by the forcing, is enforced 
according to the procedure described by Quadrio & Ricco |8j , who reported an overall error 
in determining the friction coefficient below 1%. We have also verified that the results are 
not affected by the chosen value of L x . One simulation with A+ = 1875 and large DR 
has been repeated by doubling L x at L x = 12irh (and of course doubling the number of 
streamwise Fourier modes). The value of Cj measured in the simulation with the longer 
domain differs from the one computed by the simulation with standard domain length by 
much less than the error, mentioned above, related to the finite averaging time. 

The simulations were run on a computing system available in dedicated mode at the 
Universita di Salerno, taking advantage of its large computing throughput to run several 
cases at a time. The system possesses 64 computing nodes, each of which is equipped 
with two dual-core AMD Opteron CPUs. The single computational case took about 10 
days of wall-clock time when run in parallel by using 8 nodes. Up to 8 cases can be run 
simultaneously, so that the wall-clock time for the entire study was about 7 weeks. 



IV. TURBULENT FLOW 

The effectiveness of the steady spatial forcing (T5|) in reducing the frictional drag is as- 
sessed by examining the value of the skin-friction coefficient obtained in several different 
simulations, in which the amplitude A and the wavelength A^ of the wall forcing are system- 
atically varied. The measured coefficients are then compared to the value of the reference 
(unforced) flow. The savings in driving power will then be compared to the energetic cost 
of the wall velocity distribution, in order to asses the net power saving made possible by the 
spatial forcing. 
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A. Pumping power saved 



Our simulations are performed at a fixed flow rate, so that a decrease of the frictional 
drag translates into a proportional decrease of the mean streamwise pressure gradient and of 
the power required to drive the flow. The power to drive the flow against viscous resistance 
is defined as: 



UbL x L z f 



V dr = — / (t x/ + T XtU ) dt, 

where r x ^ and r x>u are the space-averaged value of the streamwise component of the wall 
shear stress, evaluated at the lower and upper wall respectively, and (tf — U) is the time 
interval over which the time averaged is carried out, after discarding initial transients. 

The percentage saved power P sav is expressed as percentage of the power Vdr,o required 
to drive the flow in the reference case, and is defined as follows: 

P — 1 nn ^'' ~ ^ >dr 

r sav J-UU 

rdr,0 

The quantity P sav exactly corresponds to the percentage of friction drag reduction, and 
is expected to be a function of A and A^. 



In figfj P sav is p 
A + . Available data 



otted first as a function of A+, for different values of the amplitude 
obtained with temporal forcing are also plotted, with the oscillation 
period translated into an oscillation wavelength through the near- wall value of the convection 
velocity, namely = 10. 

The capital information drawn from this plot is the validation of our working hypothe- 
sis: the present steady forcing indeed parallels the unsteady oscillating-wall technique when 
oscillation period and forcing wavelength are related through U w . In analogy to the oscil- 
lating wall, that yields the maximum DR at a well-defined oscillation period T + , namely 
Topt — 100 — 125, the steady forcing yields the maximum DR at a well-defined wavelength 
A+, namely A+ opt = 1000 — 1250. It is striking how well the prediction X x opt = U^T^ pt 
is confirmed by our simulations. DR is observed over a very wide range of wavelenghts, 
200 < A+ < 8000, analogously to the oscillating wall. 

The effects of the SSL on the turbulent flow are, however, only qualitatively similar to 
those of the TSL. Quantitative differences do exist, and in particular at a given wavelength 
the spatial forcing is observed to attain larger values of DR when compared to the temporal 
forcing at the same amplitude and equivalent period. FigJS] shows P sav as a function of the 
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FIG. 4: Percentage drag reduction P sav as a function of the forcing wavelength. Empty symbols 
are temporal cases from Quadrio & Ricco [g]: for them A+ = U+T + (see text). 

amplitude A + , for different values of A+. The plot contains data from spatial as well as from 
temporal forcing cases. The dependence of DR on the forcing amplitude is qualitatively 
very similar for TSL and SSL: DR monotonically grows with A, but the increase saturates 
to an apparently asymptotic behavior. Again, the space-dependent forcing appears to be 
more effective in reducing turbulent friction: for a given amplitude the spatial forcing yields 
a 20 — 30% larger DR than the temporal forcing. A maximum DR of about 52% is observed 
at A+ = 1250 with a forcing amplitude of A + = 20. 
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FIG. 5: Percentage drag reduction P sav as a function of forcing amplitude. As in figJU empty 



for comparison recall that T + = Aj/Wj 



symbols are temporal cases from Quadrio & Ricco 
(see text). 

B. Power expended at the wall 



In addition to Vd r , the forced channel flow has an additional power input V req , that is 
required to enforce the wall motion against the spanwise shear stress. V req is defined as: 



req 



W (t z p + 7V„) dxdt, 



(9) 



tf - U Ju jo 

where t z is the space-averaged value of the spanwise component of the wall shear stress, and 
W is the spanwise velocity of the walls. 

The required percentage power P req is defined in terms of the friction power Vd r ,o of the 
reference flow as P req = V req /Vd r ,o- Preq is presented in figj6] as a function of the forcing 
wavelength. Again, spatial as well as temporal forcing data are included. Of course P req 
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FIG. 6: Percentage required function of the forcing wavelength. As in figlH empty 

symbols are temporal cases from Quadrio & Ricco [8]. The line-connected asterisks represent P req 
obtained from the laminar solution with A + = 12, from formula ([5]) for the TSL and from formula 
© for the SSL. 



assumes negative values, i.e. work has to be done against the fluid viscosity. Comparing the 
two forcing methods, one can easily appreciate how the spatial forcing presents an energetic 
cost that is smaller than the cost of temporal forcing, by approximately a factor of 2. 

The line-connected points in fig J6] represent the values P req that can be computed from 
laminar theory. TSL and SSL are considered at A + = 12, for which formulas (jSJ) and (J7J) 
respectively are analytically integrated. The good agreement was expected, since it was 
already observed (see figj3]) that laminar and turbulent mean profiles of spanwise velocity 
are coincident. 

Lastly, a net percentage power saving P net is easily defined by comparing P req and P sav , 
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FIG. 7: Percentage net power saving Pnet a function of the forcing wavelength, 
as follows: 

Pnet Psav Preq- 

The net power saving P net is plotted in figJS Since it has already been observed how SSL 
yields higher DR at lower energetic cost, of course here a significant difference in net gain is 
expected when comparing SSL and TSL. P net shows indeed an interesting maximum value of 
about 23% at A + = 6, while remaining positive over a wide range of wavelengths. Moreover, 
positive P net are found for rather large amplitudes: at A + = 12 a net gain of about 5% can 
still be observed, whereas the TSL at the same amplitude presents a net loss of about 30%. 
It is worth noting that the search for the maximum of P ne t cannot be considered exhaustive, 
and thus the presently observed maximum value of 23% at A + = 6 should be regarded as a 
starting point for a refined search. 
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FIG. 8: Isosurfaces for the quantity A^ = —0.03. Flow is from left to right. Top: reference flow; 
bottom: spatial forcing with Aj" = 1250 and A + = 12. 

C. Flow statistics 

Observing the main statistics of the flow and even a few instantaneous snapshots may help 
understan ding how SSL affects turbulence. (Additional statistical quantities are presented 
elsewhere 124)- In figE] isosurfaces are visualized for the A2 quantity introduced by Jeong 
& Hussain [25[ and since then often used to identify turbulent vortical structures. The level 
is set at A^ = —0.03. The top plot is for the reference flow, and the bottom plot is for the 
same flow subject to the effects of the spatial forcing, with A + = 1250 and A + = 12. This 
is one of the cases with the highest drag reduction, i.e. 45%. The SSL clearly appears to 



16 




FIG. 9: Streamwise mean velocity profile in law-of-the-wall form. Dotted line is the linear velocity 
profile («)+ t = y+. 

modify the near-wall turbulence and its structures. 

Similar considerations can be drawn from single realizations of the flow as well from its 
statistical description. Figj9] reports the mean velocity profile in the law-of-the-wall form. 
The modification to the profile for the TSL and SSL cases are analogous when compared 
to the reference flow, but the effects are larger for the SSL. The drag reduction manifests 
itself through the thickening of the viscous sublayer, that results in the upward shift of the 
logarithmic portion of the velocity profile, as previously documented for other DR techniques, 
for example riblets j^ . 

Another relevant statistical quantity that is modified by the action of the SSL is the 
turbulence intensity. Fig{TU] presents the r.m.s. value of the fluctuations of the streamwise 
velocity component. To examine the effect of the lower Reynolds number implied by the 
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FIG. 10: Root-mean-square value of the fluctuations for the streamwise velocity component. 

strong drag reduction in a constant-flow-rate simulation, this plot is still scaled using inner 
variables, but both the friction velocity of the reference flow and the actual friction velocity 
of the drag-modified flow are used as velocity scale. The strong reduction of fluctuations 
is deemphasized with the latter form of scaling, and what remains is the structural effect 
of the SSL after subtracting the effect due to a smaller Re T . The curves however remain 
significantly different; the main effect is a decrease of the fluctuation intensity, together with 
a displacement of the position of the maximum intensity further form the wall. Similar effects 
have been already documented for other DR methods. In particular the same observation 
has been put forward 8J for the oscillating-wall. It may be useful to remind the reader that 
effective DR methods exist, for example the active opposition control [27|, where u^ ms is 
unaffected when properly scaled; on the other hand, as discussed for example by Jimenez 



281 ] . experimental evidence exists that wall roughness may reduce the near- wall peak of 
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turbulence intensities while increasing drag. 



V. DISCUSSION AND CONCLUSIONS 

This paper has studied a new form of boundary forcing for wall-bounded turbulent flows, 
that consists in imposing at the wall a steady distribution of spanwise velocity, modulated 
in the streamwise direction. In this study only sinusoidal modulation has been considered. 
Main motivation was to find a steady counterpart to the oscillating-wall technique. The link 
between the two kinds of forcing is the convection velocity of the turbulence fluctuations, 
that takes a well-defined non-zero value U w at the wall and is thus capable of transforming 
a time scale into a length scale and viceversa. 

Thanks to a number of direct numerical simulations, the behavior of this new forcing in 
the parameter space has been determined, and DR up to 52% has been observed for A + = 20 
and = 1250. For all amplitudes, the forcing wavelength that yields the maximum DR has 
been found to correspond to the optimal period of the oscillating wall converted in length 
through U w , thus confirming the validity of the analogy between temporal and spatial forcing. 

This analogy has been extended further by studying the laminar case: this was known to 
be relevant to the oscillating-wall technique, since in the turbulent case the spanwise profile 
after space-time averaging is identical to the laminar solution. The laminar solution, that 
can be written in terms of Airy functions, has been determined for the spatial case too, 
and it has been additionally verified that the turbulent spanwise flow when phase-averaged 
is identical to the laminar solution. This property can be leveraged to predict the power 
required for the control. We hope that a predictive quantity, similar to the parameter 
discussed by Quadrio & Ricco js, 21] and capable of describing the DR effects of the forcing, 
could be envisaged on the basis of this analytical solution. 

Together with qualitative analogies, there are quantitative differences between temporal 
and spatial forcing. The spatial forcing is more efficient in terms of DR, from the point of 
view of both absolute DR and net power saving. In particular, a net power saving as high 
as 23% has been computed at A + = 6, with one unit of forcing power translating into 3 
units of saved power. This is more than 3 times the largest net saving documented for the 
oscillating wall, and also significantly larger than the benefit obtained with passive devices 
like riblets, that are reportedly capable of a saving up to 8-10% in laboratory conditions 
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29]. 

The present form of spatial forcing is certainly realizable in principle, and an experi- 
mental setup is indeed under construction that will help shedding light on effects like the 
dependence of DR on Re. At the same time, we do not consider the present forcing, though 
steady, directly suitable as yet for practical applications. However, the successful design of 
a steady control law is one important step towards the realization of a passive drag-reducing 
device. In this framework, the sensitivity of the turbulence near-wall cycle to a well-defined 
streamwise lengthscale is a fundamental result, that paves the way to the search for an effi- 
cient drag- reducing type of roughness. In this context, we have recently applied for a patent 



301 ] concerning a design method of roughness distributions yielding improved benefits over 



conventional riblets. Such distributions may contain riblets with sinusoidally ondulated 
crests with wavelength X opt . A later paper by Peet et al. 31] has confirmed by numerical 
simulations the potential success of similar sinusoidal riblets, reporting a 50% improvement 



over conventional riblets. Since in 



31] the wavelength was in the optimum range, but only 



a single amplitude was tested, we are confident that even larger benefits can be attained by 
such a purely passive technique. 
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VII. APPENDIX 



In this Appendix an analytical, approximate solution of the Navier-Stokes equations for 
the laminar flow between indefinite plane walls is described, where the non-homogeneous 
boundary condition 

w(x, 0, z,t) = A cos(kx) (10) 

is imposed to the spanwise velocity component. The boundary forcing creates a layer of 
alternate spanwise motion, which develops close to the wall and resembles the temporal 



20 



Stokes layer. We have referred to it in this paper as the Spatial Stokes Layer (SSL), in 
comparison to the conventional Temporal Stokes Layer (TSL). An analytical solution will 
be derived now for the velocity profile of the SSL, under the assumption that the wall-normal 
length scale characteristic of the layer is small in comparison to the channel width. 

The solution is steady, and thanks to the spanwise invariance of the differential system, 
including its boundary conditions, all the z derivatives can be dropped from the momentum 
equations, which then read: 



du du 1 dp ( d 2 u d 2 u \ 

u r x +v Yy = -p^ + v \d? + W l (lla) 



dv dv 1 dp ( d 2 v d 2 v s . 

u di + v Yy = --plFy +v {a? + d? ] ' (Ub) 



UUU UUU I V UU V uu \ , . 

u ~ + v — = z/ (lTT + irT ) • ( llc ) 



dw dw ( d 2 w d 2 w" 
dx dy \dx 2 dy 2 
This highlights that Eqs. flllal ) and ( lllbl) decouple from Eq. (lllc| to form an independent 
two-dimensional problem, unaffected by the inhomogeneous boundary condition ( 1T0]) . Its 
solution is thus the classical laminar Poiseuille solution, that gives a parabolic longitudinal 
velocity profile, and predicts a wall- normal velocity v = everywhere. Thus Eq. (Illcj) can 
be further simplified as: 



u w u uu u uu , . 

U ~ = + U TTY, (12) 



dw d 2 w d 2 w 

— = V 7T + V 5 

dx dx dy 

In ([12]) u = u(y) is the (known) parabolic Poiseuille profile. The PDE (fl2l) is linear, and 
thus in the following we will consider A = 1 without loss of generality. 

At this point a boundary-layer approximation is introduced. We suppose the SSL to be 
confined in a thin region close to the wall, and to vanish at a distance definitely smaller 
than h. If 5 X indicates the characteristic thickness of the SSL, we are requiring that 5 X <ti h. 
Even before giving a precise definition of 5 X , we have seen in fig f5] that this requirement is 
satisfied when the flow parameters are set within the range of interest (which is where the 
spatial forcing achieves a substantial DR in the turbulent regime). 

For small S x / h, u(y) in Eq. fll2|) can be replaced with the first term in its Taylor expansion: 

u(y) w u yfi y. 

If A x is comparable or larger than h, the boundary- layer approximation implies also that 
d 2 w/dx 2 in Eq. (fT2"]) is negligible compared to d 2 w/dy 2 . We are thus left with the well-posed 
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problem: 

with boundary conditions: 



dw d 2 w . . 



w(x,0) = cos(kx) (14) 
lim w(x, y) = 

y^oo 

to be solved in the domain y e (0, +oo). 
Its general solution w(x,y) has the form: 

w(x,y) = M[e i ™F(y)], (15) 

where the function F is complex valued, F(y) : R — > C. 

By substituting the functional form ([15]) into Eq.( |T3l) an ordinary differential equation 
for the unknown function F is obtained: 

iv- X Ku y »yF{y) = (16) 

Its boundary conditions follow directly from (fill) : 



»[F(0)] = 1, 9f[F(0)] = 0, lim F(y) = 0. (17) 

To simplify notation, the factors v~ X Ku y $ in the l.h.s. of Eq.t JTBT) are written in terms of 
a single parameter S x , which has dimensions of a length: 

I/ _1 K« V)0 = 5" 3 , (18) 

Eq. (fTB"]) then becomes: 

d 2 F 

*?V F = ( 19 ) 

Introducing the change of variable y = i5 x y, and redefining the unknown function as 
F{i5 x y) = F(y) turns Eq. (|T9|) into the following Airy equation: 

d 2 F 

mv) = w - (20) 

Infinite solutions of an Airy equation exist for y spanning the whole complex plane, when 
derivatives are considered in the sense of analytic functions. These solutions are linear 
combinations of the two special functions Ai(y) and Bi(y) 

F{y) = aAi(y) + PBi{y), (21) 
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which are called Airy functions of the first and second kind respectively. The general solution 
fl2Tj) has an alternate representation |20|, which turns out to be useful in our case: 



F(y) =~fAi(ujy) + 9Ai{uj 2 y), 



(22) 

where uo = e~ i2//37r . It can be shown [20( that among Ai(y), Bi(y), Ai(uy) and Ai(u 2 y), the 
only base function satisfying conditions (1T7|) (up to a normalization factor) is Ai{uo 2 y). As 
a consequence, the solution that satisfies the required boundary conditions is: 

Ai(u 2 y) 



F(y) 



Ai(0) 



(23) 



By substituting this solution into (|T5l) . the expression for the unknown function w(x,y) 
is eventually derived: 



w(x,y) = CJR 



(24) 



where C x = Ai(O) -1 . 
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